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THERMAL MODE! OF LASER-INDUCED SKIN DAMAGE: 
COMPUTER PROGRAM OPERATORS MANUAL 


INTRODUCTION 


The Skin Medel is a mathematical model that predicts the 
extent of thermal damape and degree of burr produced in skin 
by laser radiation. The model was developed by IIT Research 
Institute and is a result of many years of improvements ia 
thermal damage modeling techniques. The mathematical basis 
for temperature predicticns is the following heat-conduction 
equation in cylindrical coordinates: 


a) Pas 
Coe = q(r.z.t) + 


Kooy 5) 34 e) av 
rove * Te (RRe) + 5g (gD . 


nd 


- B(z)C, - q(r.z.t) 


= density of tissues, g/em? 

= specific heat of tissues, cal/p-°C 

= blood flow, g/em-sec 

= specific heat of blood, cal/g-°C 

thermal conductivity, cal/cm-sec-°C 

= rate of heat deposition from laser. cal/em>-sec 


= rate of heat loss in transforming water into 
steam, cal/cmj-sec 


where 


i2] 
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= radial distance, cm 
= time, sec 
= temperature, °C 


nu cana 


= axicl distance, cm 
In addition, equations presented on the next page provide for: 


e heat losses to the ambient air, 
e steam blisters, and a 
e hair follicle. 


Initially, the skin is considered to be at a uniform 
temperature measured with respect to the surrounding air. 
Biological heating is neglected. Irradiation from the laser 
is reduced as a consequence of interceptions of radiation by 
hairs and reflection eee the surface of the skin. In this 
regard internal reflections ware considered to occur at the 
surface. 
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Onee the radiacteot enters the skin ita intensity, q, is 
Moen ted qeeordin: fo Reer'’s Law, namely 


a en oe action f ea) (2) 


Vo 
what ou radiant oatensi py considered a: ao step function 
in time, cal/em -see 
\ absorption coefficient considered as a step 


function with respect to depth «2, L/em. 


Beam diverpence is neglected. Equation 2 is then used to 
determine the rates of energy deposition q(r,z,t) for equa- 
tion l 


Blood flow B, specific heat C, density “, and thermal 
conductivity K are considered to vary with respect to depth. 
Blood flows remain constant with respect to time and act aa 
a heat sink. Blood is not treated as a means of heat trans- 
for Thermal properties and density vary with the amount of 
water in the tissues 


When water commences to boil, the resultant water vapor 
is allowed to escape from a given depth of skin. Resultant 
rates of heat loss 7 per unit volume are then introduced into 


equation 1. Additional heat losses are allowed to occur from 


the surface of the skin by natural convection/radiation. 
These losses are described by 


av 
-K=— | = Kv} 
dz i 3 
z= mea 2) 


where H = heat-transfer coefficient, cal/om*-sec-*C. 


This coefficient is assipned values according to whether or 
not sweat is initially present. The coefficient is not 
chanred with cime. 


Following water evaporation, densities, specific heats, 
and thermal conductivities are decreased according to the 
amount of water lost at specific locations. In addition a 
steam blister is allowed to form at the base of the epider- 
mis provided it reaches the temperature DTEMP. Once a blis- 
ter forms, thermal resistance is introduced to account for 
the presence of water vapor. 


Mathematically this is given by 


dv 7 ; 
aS Pe. sce: H(v| vi ). (4) 
hoe! z +s z -S z +8 
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Tha aa atdem omen deat the haat capactty ot che water vapor 
ae that heat fluses crosaite the two blister nurfaces are 
vader teal These tluaes ate expressed by the product of che 
Peat ttanaster caettiotent dH. acrom the bltater and the tem 
woriture dft( ference between the sur facer Deprh of blister 
formation fs desertibed be oy. and representa a vaniahing iv 
amall dfspiacement of the upper and lower blister surfaces 


The model provides va hair tebldele atoaopiven depth on 
the ands ef the Laser beam Radiation is considered to be 
entitely absorbed by the folltete Meat deponition rates 
direct ve below the ftolblttiele are sere 


The heat -conduct ton equation is ee by finite 
differences and then solved with an explicit- implicit 
alternating-direction technique developed bv DW Peaceman 

and WH. UW. Rachford (1). This technique solvea tie finite- 
difference equations explicitiv in s and tmplicitty in r for 
odd time steps, and implicitly in s and explicitly in or for 
even time steps. In explicit: calculations, existing Cempera- 
Lures are used toa represent thermal pradients, in daplicit 
calculations, future temperatures are used This approach 
reaults in a set of equations that are solved uaing ordinary 
matrix alpebraoo Larger time intervals can be used with thin 
technique than with standard explicit: finite-difference methods 


The model ages the predicted temperature risen to deter- 
mine irreversible tlasuc damare os applying Henrique'’s damape 
equation fer temperatures above a yiven value 


Vodrye yt) CpexplC/T (rc tyyat (5) 


where \i(r,2.0) © incremental damage at point », + 
C, and C, = censtantes riven bv arrav DAM in nomenclature 


TGQ = absolute temperature 
At + inerement of time. 


Irreversible tissue damiye ie defined aa vecurring whenever 
the intesral of Ac over time iso sreater than or equal te 

1. From this mathematical basis, the model has the capability 
of predicting temperature rises, damage thresholds, the ex- 
tent of damape within the skin, and degree of burn based upon 
the resultant ¢ values and peak Cemperatures. 


1. Peaceman, Do OW. and Ho OH. Raehford. Jr The numerical 
solution of parabolic and elliptic ditferent tal 
equattions. 1 See Indust Appl Math & 78-41 (1955) 
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(estened Poa manimup tlenibility. tin model offers wide 
vatiabilitv in characterizing tissues and laser exposures. 
it accommodates variations in laser radiation characteriatics 
awd in optical. thermal. and physiological propert tes of the 
abba The medel” oo desi: enables the user te print out analy 
these partians of the output information which he destice. 


The purpess cot this manual (ds to give the user a basic 
wnderatanding ef the wdel'sa capabitithes and how to uae it 
within the Limits of those ¢ ilictes. Additional ir:forma- 
tlen en the code can be obtained from the LITRI Technical 
Report “laser-induced Thormal Damage of Shin" (7). 


The manual brietiv deacribes (1) the .apabilities atid 
bimitations ef tho model an they pertain to the source, the 
shin, the mechanics of the program, and the output deaired, 
and (27) the basic input required, listing the required cards. 
their order. and app opriate formate and optians available 
rv the user in regard to printed aut put Appendix A presents 
a ploasary of all input parameters. 


CAPABILITIES AND LIMITATIONS 


The user is responsible for adequately describing the 
exposure condicions to be aadeled. must describe the 
incident Su aaa the media, the mechanics (temporal and 
apatial grid) of the program. and the output dealred. This 
section presenta a broad overview of che capabilities and 


Limitatdona aa they pertain to these areas. 


In develbopine the propram, several major asaumptiona are 
made First. the skin veonetrices are almulated in cvlindrical 
coordinates Second, all cetleetions occur at the surface of 
the skin Third, beam profilee are identical at all depths 
within che skin, and are radially svametric Fourth, therwal 
Properties change ontv with water Logs and not with ¢ empera- 
ture Fitth, Slood flows act solely as a heat sink. Sixth, 
steam blisters remain intact. while water beneath the blisters 
becomes a superheated liquid without deforming tissues. 


The model haa oc number of features which give the user 
flextbilitv in describing che incident radiation jn terme of 
ite apatial, spectral, and temporal properties. The spatial 
profile of the beam may be designated as uniform, gaussian, 
Or irregular. For untftorm or gausaian profiles, the user 
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apecities the bewn tafius and the incident power For ir- 
rervular prefiles, the user constructs the desired beam pro- 
files by apecifvins the tneident power and the relative 


intensities snd smsociated radial distances fram the center 
et the been 


temporal varfations of the incident radiation are speci- 
fted by selecting the duration of the pulses, the repetition 
rates und intensitv, and the number of pulses. The aodel 
assumes all pulses te be square with respect te time How- 
ever. temporally varving pulses can be constructed from a 
collection of back-to-back square-shaped pulses. The pulses 
awe ‘spe rates nay be considered identical (noncoded) 
or variable (coded) In that the model ogly predicts thermal 
damaye. pulae duratioas should exceed 10°:' seconds. For 
shorter pulse durations, other damage mechanisas may predom- 
inate 


Skit media are described in ternss of optical, thermal. 
and phyasfologtcal characteristics. The optical parameters 
include coeffietents for absorption and reflection. Thermal 

araaeters include heat capacity. thermal conductivity, and 
load flow associated with individual skin and subcutaneous 
lavers plus their inittal temperature. The skin is considered 
to be composed of two layers, namely the epidermis and dermis. 


A single homogeneous laver is used to represent subcut eneous 
tissues. 


tially, the thermal jroperties and density are con- 
stant within che dermis and within the subcutaneous tissues 
and varv linearly acroas the epidermis. paca 44 rties 
and blood flow rates are constant across layers ich are 
not necessarily coincident with the epidermis or dermis. 


The model automatically selects the time steps. However, 
the user specifies the aize of the increments for the spatial 
grid in the skin media. The dimensions of the apatial and 
temporal tnerements represent the limits of resolution of 
the mode! predictions. The spatial grid has unifcrm incre- 
ments in the region of highest temperature rise and constantly 
expandin, increments .way from the highest temperature regions 
The temporal yrid has constantly expandiny time steps from 
the berinning and end of each pulse. 


The user selects the range of spatial coordinates used 
to print the temperature riscs and damage. ilso, the user has 
the opcion to print and plot cemperature rises at any selucted 
time. A separate plotting routine uses input data cards 
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desired Ali porametets required by the user as fuput are 
detined in Arpendix A Program Tisting te presented tn Appen- 
fix Ro When selvidns a problem, the user must describe the 
iteetsemt radiation, and must also specify the parameters 
voverming the mechanical operation of the program, such as 

thw kice of the grfd required and the output deatired. 


Te simulate the radiation incident on the skin media, 
the user tas the option of specifving a uniform. yaussians, 
or atrentiarsbeam irradiance profile.  Axdal symmetry for 
all beam profiles da assumed For uniform-beam profiles. 
the beam radius CRUNTFY must be specified. For saussatan-beam 
prefiles, RIM must be specified at a particular relative ine 
tensity pedag cc For Lrrepular-beam profiles, the abso- 
lute or relative drradiance profile must be specifted on a 
peint-by-point basie by listing che irradiance value, PX(L), 
at the radial distance, RXCLY, from the center of che beam. 
The tetal number of specified irradiance points is equal to 
LR The model will interrate the beam profile at radial ine 
tervals CRINTY. to arrive at a series of conatant intensities 
centered about each radial prid painr. 


The model can treat coded pulses which differ in inten- 
sity and durat don or nonceded pulses which are identical as 
are Che times between pulses Coded pulses are considered 
te fellow one another with no clapsed time between pulses, 
Single pulaes are considered as noncoded. Data cards required 
for pre; ramming are described in Table 1. Further informa- 
tien is presented ar the bottom of Table 1. 


Inport cards are read in sequence according to the card 
numbers presented in the Listing and Table 1. Tiis table 
summardces the format statements along with che input parime- 
OUT, 


ALY cards cited in Table lb. except for card number 21 
Ath these toblowed be terters. are needed for each computer 
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| 
TABLE DL ENET? DECK i 
: ae t Ne Format Parameters 
bs UP, 10R? } LBL ABSW.COML.COMZ SHO. SHI. SH2,WO. Wl 62,0) | 
F e 10E? ) DO,D1,0).05.28L. WW. OTEMP, DHF. SHE 
lor? TALL) Led kp 
i + 10K? 4 CH, (BCL) Lt, LB) 
s 10F? TSWEAT DEPID, TEPID. TDERM. TO. TE, ZDEP 
oe Lor? 4 DANC 2.1) DAM(1.2).DAM(2, 1), DAM(2, 2) j 
} . lOF? 4 TOL(L) Le) .27 4 
i as 17,1087 3 LZ. ABSC.H1.N2, HAIR. REF } 
| 99 10F? « TH(L) Lek. L2 1 
10 107 4 ABS(L) .L=3.L2 i 
t ¥ 1017 TPROF . MK 3 
DUMIFS OF? 8.207 RUNIF ttn | 
S.GAUSS" QF? 4.217 RIM.CUT.NL.N F 
12, URREG* = 1017 LRN LN 
12. TRREG@ = 1087.3 RX(L) Le 2. UR 
12. 1RREG® = 1087 3 PX(L) Let La i 
13 GF7.1,617 WAVEL. XC. 2.22. LASER, NPULSE ,.NGX. NTX ; 
ta.nc.h 17, 1087.3 NTP. DPULSE. POWX 
1e.nc xh 17, 1087.3 -- REPET.. DPULSE. POWX { 
l2.nc.GPe =. 1017 NG 
| la.uc.crp® = 1087 NPG(L) ,L=1,5G ; 
| le.wc.cP> 1037 WPR(L) .Le1.Nc i 
te,copep> = =—s 017 NTP j 
le.copep> = Lor7.3 DPULSC(L) ,L=1, MPULSE i 
| la.copep> =-1087.3 POUERC(L) Lol ,MPULSE | 
1$ 1017 ¥ i 
le 1017 1D2, 1D2,JD1,JD2, ITYPE. KTYPE , KX ‘ 
1? 1087.3 TIMER (L) .Lel RX 
18 10E7..3 TIME2(L) Lol AX 
19 17,1087.3 0 NW, (ZW(L) Let Ww) 
20 10E7.3 DB]. DB2,DB3. 085 
| 2 1? TUL LIZ 3.1 Js? 
| 
| 


For single-pulae exposures use card la,NC.1, for multiple 
noncoded pulses use cards 14.NC.X, for multiple noncoded 
Pulses to be grouped use card 14,NCX, ond cards 14 .NC GP 
(ave text). For coded pulses use cards 14,CODED. 


run. Card 2) is required only when one wishes to punch tem- 
perature rises on cards. The other exce,tions are cards 12.... 
(beam profile), and cards 14....(pulse torulation), For uni- 
furm profiles use card 12, UNIF} for gaussian profiles use 
card 12, GAUSS; and tor irregular profiles use card 12, IRREG. 
A description of the cards 14.... needed for describing pulse 
characteristics is presented at the bottom of Table 1. In 
this reyard, provisions have been made to minimize computer 
exccution times caused by long noncoded pulse trains. This 
provision is optional and requires use of cards 14, NC, GP. 
These cards result in the replacement of various sets of con- 
secutive pulses by a "single pulse" which extends from the 
start of a given pulse to the start of the pulse immediately 
tollowing each set of consecutive pulses. Power is set equal 
to the total energy associated with the pulses divided by the 
duration of the "single pulse." By this means total energy 
is conserved. Resultant power is termed mean power. 


To introduce mean powers into long pulse trains, one 
should first select several sets of pulses spaced at regular 
intervals along the length of the pulse train with at least 
two individual pulses between the sets of pulses and at the 
end of the train. Incremental damage is only predicted during 
individual pulses hee pole the first of each series of indi- 
vidual pulses. Weighting factors are used to account for 
damage produced during remaining pulses. 


An example of the above procedure is shown in Figure 1. 
Here we have chosen a pulse train consisting of 40 pulses. 
Number of pu’ses replaced by mean powers is indicated by 
numbers over the shaded area. Table 2 presents the values 
needed for the arrays NPG(L) and NPR(I.) to structure the 40 
pulses in the fashion illustrated by Figure 1. 


TABLE 2. PROGRAMMING 40 PULSES IN MANNER 
IL].USTRATED BY FIGURE i 
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L NPG(L) NPR (L) L NPG(L) NPR(L) 
1 2 0 12 1 ? 
2 1 0 13 5 0 
3 1 ? 14 1 0 
4 5 0 15 1 7 
5 l 0 16 3 0 
6 1 ? 17 1 0 
? 5 0 18 l 2 
8 1 0 19 1 1 
9 1 7 20 1 1 

10 5 0 21 1 1 

11 1 Q 40 40 
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The index L indicates the order in which individual 
pulses or mean powers are to be considered. The array NPG(L) 
represents both the number of pulses and how they are to be 
treated. For normal treatment of pee NPG(L) is sat equal 
to 1. NPG(L) values greater than 1 result in replacement of 
the specified nunder of pulses by the mean laser power. 


The array NPR(L) represents the number of pulses repre- 
sented by individual pulses shown as shaded in Figure ?. 
These pulses have L values of 3, 6, 9, 12, and 15 in Table 2. 
Each of these pulses represents 6 additional pulses--three 


immediately preceding it and three immediately following it. 


All other NPR values should be zero. The sum of the two 
cheers ure ane NPR should equal the total number of pulses 
n the train. 


Thicknesses of the epidermis and dermis will vary widely 
with individual and location. Epidermal thicknesses for 
humaas range from about 0.005 to about 0.015 cm while dermal 
thicknesses range from about 0.15 to 0.30 cm. Recommended 
values presented in the Glossary represent average values for 
pig skin. 


Optical property data are presented in Table 3. These 
data were acquired from several sources as indicated in ref- 
erence 2. For those coefficients not presented in the table, 
use the coefficients for water given in Table 4. These 
values are not critical in that at long wavelengths much of 
the irradiation is deposited within shallow depths of the 
epidermis for which coefficients are known. Recommended 
values for the remaining input are presented in Table 5 and 
Appendix A. 


Table 6 is presented to assist the user in identifying 


parameters requiring change in altering various features such 
as wavelength, laser profile and pulses, and skin. 
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TABLE 3. OPTICAL PROPERTIES OF IRIAN ERIW (ref. 2) 

F 

abeorption ove tfictente : 

et) 0 200 “0 22 ry) “ 

er 0 110 2 20 39 

30 sxe 10 0 $5 13 13 19 . 

rv nt) 0 29 12 12 13 4 

oe » 4 ° 25 e ua 12 4 

a ee ? 0 40 10 a2 a 

0 808??? 0 te 10 10 10 : 

0 806378 0 ) uw 10 10 : 

1.00 39 98 0 298 rT) 10 10 F 

1.10 se 8 ) an 13 10 10 EH 

1.20 38 40 1 233 12 12 12 4 

1.423932? N 217 a1 24 yy 4 

1.70001d ? 231 20 2? 25 4 

es | a3 205 rr) 28 35 : 

2.20 » 3 7 202 0 23 ay } 

2.40 » 3 30 24 50 28 54 : 

3.00 2 2 386 Yve : - : 4 

4.00 2 2 145 380 . : - 

4.? 2 2 420 $72 - - : J 

3.00 2 2 312 401 : 7 - 

f 6.00 2 2 2241 1290 - - - 4 

E 7,00 es 574 999 F - : d 

8.00 a $39 921 : : : ‘ 

; 9.00 too 337 821 - . - 4 

: 10.00 1 1 638 727 - - - 1 

11.00 1 1 1106 911 : - : F 
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TABLE 4, 


0.200 6.9115-10°° 
0.225 2.7987 107? 
0.280 1.6839 1072 
0.275 3.0739 -:10°2 
@.300 6.7021 20°> 
0.325 4.1759:20°> 
0.380 -2.3298-10° 
0.975 —-1.1729-207? 
0.400 $.6434-10°4 
0.425 3.8438-1074 
9.430 2.8484-1074 
0.475 2.4796 1074 
0.500 2.5193 -10°4 
0.525 3.1595 -1074 
0.580 4,4782-107* 
0.575 7.8676 207° 
0.600 2.282910? 
0.625 2.7948-107° 
0.650 3.1706-107° 
0.675 4.1516 -107° 
0.700 6.0139 107° 
0.725 1.$840-1072 
0.780 2.6138 -10°7 
0.775 -2.3998-1072 
0.800 = 1.9635.1072 
0.825 2.7722 1072 
0.880 4.3317 -1072 
0.875 $6154 -1072 


ABSORPTION CONSTANTS OF WATER (ref. 4) 


6.7858 1072 
1.4400 -207 


4&. Hale. G. M. and M. R. Querry 


of water in the 
region. Appl Opcics 12(3): 


200 nm to 
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0.950 3.8757 107! 
0.975 4.4852 10° 
1,000 3.6317-107! 
1.200 1.0357 
1.400 11,2387 10°" 
1.600 6.7152 
1.800 8.0285 
2.000 6.9115 -10*¢ 
2.200 1.6508 -10°! 
2.400 5.0056 :10°! 
2.600 1,5321-10* 
2.650 31772102 
2.700 8.8430°1 
2.750 2.6961 :10° 
2.800  $.1612-10 
2.850 8.1572-10° 
2.900 1.1613 -10° 
2.950 1.2694-10° 
3.000 1.1394 10° 
3.050 9,8883 -10° 
3.100  7.7830-10° 
3.150 5.3856 -10° 
3.200 3.6285 -10° 
3.250 2.3586 -10° 
3.300 1.4023 107 
3.350 9.7905 10" 


7.2072 107 
4.8080 10" 
3.3750 +20" 


3.600 1.797720 
3.700 -1.2227-2 
3.800 1.2244.2 
3.900 1.22446 +20 
6.000 1.44511 
6.100 1.72252 
4.200 2.0588 107 
4.300 2.46942 
4.400. 2.9427 -20 
4.500  3.7420-2 
4.600  4.0158-2 
4.700 &. 1977-2 
4.800: 3.9270-107 
4.900 3.51351 
5.000 3.2165 -107 
5.100 2.7350 -207 
5.200 2.4408 2 
5.300 2.3236 +207 
5.400 2. 3969-107 
$.300 2.6504-1 
5.600 3.1865 -107 
5.700 4.4754°1 
5.800 7.1498°2 
5.900 1.3248 +1 
6.000 _2.2410-10° 


. Optical constants 
200 im wavelength 


555-564 (1973). 
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AX) See Gee ae 
e100 -2,6987:10° 
6.200 1. 7830-10 
6.300 1.137010? 
6.400 8.8161-107 
6.500 7.97485 +107 
6.000 6.7782 -107 
6.700 6.3207 +107 
6.800 6.0430 :107 
6.900 $8643 107 
7.000 5$.7446-+107 
7.100 6637-107 
7.200 $6025 +10? 
7.300  $.3430-107 
7.400  §.3020-107 
7,300 $.4622 -107 
7.600 $.4234-107 
7.700 $4019 -107 
7.800 5.397110? 
7.900 5.3924 -10? 
8.000 $3878 107 
8.200 3790-107 
8.400 5.4006 -107 
8.600 $4357 -107 
8.800 5.4978 +107 
9.000 $.5722-2107 
9.200 5.6685-1 
9.400 5.7886 -107 
9.600 5.9429 10° 


TABLE 4. (Cont inued) 

Norn), ten” t 
9.800 6. 1422-10? 
10.000 6.3897 -107 
10.300, 7.9228: we 
11.000 1.1058 «0 
41.500 1.5517-10° 
12.000 .2.0839-10° 
12.300 2.6038 10° 
13.000. 2.9483 10? 
13.800 3, 1928 10° 
14,000 3.3222 -10? 
14.300 3.3626 -10° 
13.000  3.3678-10° 
18.300 3.3764 10° 
16.000 3.314410" 
16.300 3.2896 -10° 
17,000 3.2712 -20° 
17,300 3.0806 -10° 
18.000  2.9740-10° 
18.500 2.8597 10° 
19.000 2.7982 -10° 
19.300 2.6035 -10° 
20.000 =. 2.4693 107 
21,000 2.2859 «10° 
22,000 2.1306 -20 
23.000 2.0052 -10° 
24,000 1.4902 -10° 
25.000 1.7895 -10° 
26.000 1.6916 -10° 


27,000 
28.000 
29.000 
30,000 
32.000 
34.000 
36.000 
38.006 
40,000 
42.000 
44.000 
46.000 
48.000 
$0.000 
66.000 
70.000 
80.000 
9.000 
100.000 
120.000 
120,000 
130.000 
140.000 
150.000 
160.000 
180.000 
490.000 
200.300 
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1.6011 -10° 
1.5169 20° 
1.4430-20 
1.373920? 
1.2728 10° 
1.2100-10° 
1.1973 -10° 
1.1938 -10° 
1.2095 -207 
1.2237 10° 
1.2452 :10° 
1.262110" 
1.277610 
1.291819" 
1.2294 10° 
1.0340:10° 
8.5922 -107 
7.4840:107 
6.6853 -107 
6.0661 -207 
5.5083 -107 
&.9686 -107 
&.4880:107 
4.1469 -1¢ 
3.8956 +10? 
3.4837 -207 
3.3136 +107 
3.1667 -107 
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TABLE 3, TBL(1) VALUES (ref. 3) 


hy TELL), °C L 

\ 100.9 10 
2 119.6 s 
3 132 9 12 
‘ 142.9 1) 
3 152.1 le 
6 138.1 15 
? 164.2 16 
8 169.5 17 
9 174.5 18 


TAL (L) °C 


ee ee 


179.0 
183.2 
107.1 
190.7 
194.1 
167.4 
200.4 
203.4 
206.1 


L 


19 
20 
21 
22 
23 
24 
25 
26 
27 


TRL (L), °C 


208.8 
213.4 
213.9 
216.2 
218.38 
220.8 
222.9 
225.0 
227.0 


BEA ad L indicates atmospheres of pressure caused by temperature 
LOL, 


5. Handbook of Chemistry, Sandusky, Ohio; 


Handbook Publishers, 
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rates of thermal damage 
eriterta for degres of buras 

hair tollicles 

steam bilaters 

hate vower ing 

heat-tranafer coefficients 

valculation of heat loanee From cacaptng stean 
atom pressures 
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APPENDIX A 
GLOSSARY OF INPUT PARANETERS 

lai: -Absorption coefficients of undamaged skin and 
a a aneous eiceaee: Coefficients bd Pak pt vavelongth 
and are considered constant in each of layers of 
thicknesses TH(L). 
Units. L/en 
Recommended values given in Table 3 


ABSC--Absorption coefficients of irreversibly damaged skin. 
Coefficients vary with wavelength. 


Unite: iL/cm 
Recommended values given in Table 3 


ABSW--Absorption coefficients of water. Coefficients vary 
with vavelcagth. 


Units: Il/cm 
Recommended valuea given in Table 4 


BL Ct) Le LLB" Blood flow ritea cousidered constent in TB 
ayers of tissue of thicknecees TB(L). 


Unite; g/cm? “see 
Recommerded values: 81(1)=0 
BL¢2)=0,.0118 
BL(3)*0 
C3--Specific heat of bloud. 
Unita: cal/p-°C 
Recommended vaiue: 1.0 
CONL CONS 2--Constants used to compute thermal conductiviry of 
tissues as function of water content WATER(i,j). Therwal 
conductivity is given by: 


CON] + CON? WATER(1.J)/00(1.J) 
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where WATER(T..J) - grams of water/cn? of tigaue 
oot.) = density cf tissues, g/cm 


inkts csabdctae sees CO 


Recommended values CON = 1.33°1073 
COMZ = 1.16°10 


CUT--Value of normalized gauesian profile at radius RIM. 
Unite: witless 

DBL, bB2,DR3--Cumulative thermal damage at base of epiderais 
required to produce firat, second and third degree burns 
respectively. 
Units: unitless 


Recommended values: DdOBl = 0.) 


21 
DB} = 10,000 
DBS--Temperature of skin required for fifth-degree burns. 


Units: °C 
Recommended value: 400 


DAM(1.1) .DAM(1, 2) 
DAM(2 1) ,DAM(2, 2)--Constants used to compute rates of thermal 
mage of tissues. Rates of damage are functions of the 
absolute temperature v as described below: 


exp (DAM(1,1)-DAM(L,2)/v), 317°K = vs 323°K 


exp (DAM(2,1)-DAM(2.2)/¥), v » 323°K 


(mite: DAM(L,1) and DAM(2, 1), In(V/sec) 
DAM(1,2) and DAM(2,2), °k 


Recommended values: DAM(1,1) = 149 
DAM 


DEPID--Thicknessa of “dead" epidermal layer. 
Units: cm 


Recommended value: 0.0060 (half epiderma! thickness) 
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DHF--Depth of hatr fullicle If sec to unrealiatically Large 
value such aa 1,000 em, code will diactegard follicle. 


Units om 
Revconended value. 0.715 om 
OPULSC (L) Lo 1 NPULSE--Duratiuna of each NPULSE coded pulses. 


aea considered back cto back and may have any desired 
pewer. Power considered conatant during each pulse. 


Units: ‘gec 


DPULSE--Duration of noncoded pulses over which time power {se 
considered constant. 


Units: sec 


Unite: °C 


Recommended value: 131 


DO, DL, D2 .D3--Densities of sweat, oucer epidermal layer, dermis 
and subcutaneous tissues, respectively. 


Units: g/cm? 
Recommended values: DO = 1.05 
0.50 
1.03 
0.98 


HAIR--Fraction of skin surface shielded by hairs. 


4 
eacgus 


Imitsa: unicless 
Recommended values: 0 to 0.1 


HW--Heat-transfer coefficient associated with heat fluxes 
nucrose steam blisters. 


Units: cal/em-sec-°C 
Recommended value: 6.107° 


HL, W?--Heat-transfer coefficients associated with heat losses 
from dry and wet surface, respectively. 


he. PSN tay eee i i em a i 


DIEMP--Temperature at base of epidermis causing steam blisters. 
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Recommended values: 2.0°107¢ 
7,0°10 


EDL, UD! sma.test and largest I values at which temperature 
tlave and damage will be printed. 


Hl = 
H2 = 


alba! ee a awe vl 


Unica: unitless 


SEL, 212, 123-2 indices needed for punching cards for 2 or 3 D 
temperature illustrations, I2%2 and [12 represent smallest 
and lergzeat I indices at which temperatures will be 


punched reapectively. I13 represents curve to be marked 
in illustrations. 


Units: unitless 


sen aatbaailsl bea, tea adhe tela bat nt Sl Hig os 


JPROF--Index indicating shape of beam profiles. Beam profiles 
assumed to be radially symmetric. 


Unic: unitless 


Input value: 9% (uniform profile) 
1 apausesen profile) 
2 (irregular profile) 


ITYPE--Index controlling princing of temperature rises. An 

ITYPE value of 1 wilt allow temperatures to be printed 
following each time atep: an ITYPE value of 2 for every 
other time scep, etc. Temperature rises will always be 
printed following the last time step regardless of the 
value aestaned ITYPE. Additional options for priacing 
are available via the arrays TIME] and TIME2. 


metallia il is aid aa a ts! Aa dada ds 


Unit: unitless 


JD1,JD2--Range of J indices over which temperature rises and 


amage will be printed. JD2 poreenent? smallest value 
while JD2 represents largest value. 


Unic: unitless 


JJ1 ,JJ2--Range of J indices over which tempscature rises will 
be punched on cards for vrepuring temperature illustra- 


tions. 


Unit: unitless 
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KIVPE--index indicating whether or not to punch carda for 
“temperature-rise illustrations. <A value of sero will 
skip cards while a value of i will pick cards. 


Unit unitless 


Input value: 0 (to skip punching carde) 
1 (te punch cards) 


KX--Number of time intervals TIMEL(L) to Lede ey over which 
printing of temperature rises ana damage is allowed. 


Unit: unitlesa 

Recommended valuen see TIME](L) .TIME2(L) 
LASER--Parameter describing whether pulses are noncoded, or 

“coded. Noncoded pulses arc identical while coded culees 

are not. 

Unit: unitlers 


Input value: 1 (noncoded pulses) 
2 (coded pulses) 


~-Mumber of tissue layers of thicknesses TB(L) used to 
describe blood flow BLL). 


Unit: unitless 
Recommended value: 3; see BL(L) ,.TB(L) 


-eNumber of radii RX(L) used to describe relative intensi- 
ties PX(L) of irregulor laser beams. 


Unit: unitless 


Lf--Number of tissue layers of thicknesses TH(L) used to 
describe absorption coefficients ABS(L). 


Unit: unitless 

Recommended value: 13, see ABS(L) , TH(L) 
M--Total number of axial increments. 

Unit: unitless 


Recommended vaiues: 15 to 20 
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MK--Dimension of array DTX(K), POWER(K), or XP(K). Introduced 
to. prevent overflows of storage allowed in dimensioning 
arrays with time indices, i.e. DTX, POWER and XP. 

Unit: unitless 

N--Number of radial increments. 
Unit: unitless 
Recommended values: 15 to 20 

NG--Number of groups of pulses, used only with noncoded pulse 
trains. And then only when selected group of pulses are 
represented by mean power. 

Unit: unitless 

NGX--Index indicating whether or not noncoded pulse train is 

to be treated on a pulse-by-pulse basis, or with mean 


powers representing specified sets of pulses. Index 
applies only to noncoded pulses. 


Unit: umnitless 
Input value: QO (if pulses are to be treated individually) 
1 (if some pulses are to be replaced by the 
mean power) 


NPG(L) , L=1,NG--Number of pulses contained within each of NG 
sets of pulses represented by mean power. 


Unit: unitless 
Input values: see Table 2 


NPR GL) L=]_,NG--Number of pulses represented by particular in- 
v 


{vidual pulses. 
Unit: unitless 
Input values: see Table 2 
NPULSE--Total number of pulses in pulse train. 


Unit: unitless 
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NTP--Minimum number of time steps assigne:! to each pulse 
pertod. 


Unit: unitless 


Input value: ~10 for single pulses 
2 or more for multiple pulses 


NTX--Number of time steps assigned after last pulse. Should 
be large enough so no further damage occurs. 


Unit: unitless 
Recommended value: 15 to 20 


NW--Number of trials used to approximate rates of heat loss 
for transforming water into steam. 


Unit: unitless 
Recommended value: 4; see ZW(L) 
Nl--Number of uniform radial increments. 
Unit: unitless 
Recommended value: approximately one-third of N 
POWERC (L) ,L=1,NPULSE--Laser power of successive coded pulses. 
Unit: watts/pulse 
POWX--Laser power of noncoded pulses. 
Unit: watts/pulse 
ERE eek: eeR LE UUe® of irregular beam profile at radit 
>. Values need only be relative. Intermediate values 


are linearly interpolated between specified points PX(L). 


Unit: unisless 


REF--Surface reflection at given wavelength WAVEL. 
Unic: unitless 


Recommended values: see Table 1 
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REPET--Pulse repetition rate. q 
Unit: Ll/sec : 


RIM--Radius at which normalized gaussian profile equals CUT. 
Normalized profiles have an intensity of l at the axis 
of the incident beam. 


bei oe Bete le af Alt 


Unit: em 


RUNIF--Radius of uniform profiles. 
Unit: em 


RX(L\, L21,LR--Radii at which magnitude of irregular profiles 
“equals PX(L) 


(> . 
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Unit: cm 


SHF--Radius of hair follicle. 


Hair follicle located on axis 
at depth DHF. 


Radiant energy totally absorbed by follicle. 
Unit: cm 


Recommended value: 0.02 


SHO--Specific heat of non-water constituents in sweat. 
Units: cal/g-°C 


Recommended value: 0,22 


SH1,SH2--Constants used to compute specific heat of tissues 


~~ as function of water content WATER(I,J). Specific heat 
given by: 


1 Lyman Li i ite ian It te id dll at setivnaleblls Jak Misa. tke 
dha 


SH1 + SH2°Water(1I,J)/DD(I,J) 


where WATER(I,J) = grams of water/em> of ti 


gsue 
DD(I,J) = density of tissues, g/cm 


Units: cal/g-°C 


r 
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Recommended values: SHI] = 0.37 
SH2 = 0.63 


BE EE eT SO Re PEP eS” ee rn 


SB CE) ie Le Th teknesses of LB layers of tissues within 


¢ lood flows are considered constant. Blood flows 


given by BL(L). 
Unit: em 
Recommended values: TB(1) = 0.0137 
TB(2) = 0.0219 
TB(3) = 3.000 (arbitrarily large) 


TBL(L) ,L“1,27--Temperatures at which pressures of confined 
superheated water equal L atmospheres. 


Unit: °C 

Recommended values: see Table 5 
TDERM--Thickness of dermal layer. 

Unit: cm 

Recommended value: 0.1779 
TE--Temperature of environment. 

Unit: °C 

Recommended value: 22 
TEPID--Thickness of epidermis. 

Unit: cm 


Recommended value: 0.0121 


THCL) LoL LZ -Thicknesses of LZ layers of tissues within 
which absorption coefficients are considered constant. 


Coefficients given by ABS(L). 

Unit: em 

Recommended value: TH(1) = 0.0030 
TH(2) = 0.1870 


TH(3) = 3.000 (arbitrarily large 
value) 
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TIMED CL) TIME (L) i Let KX- “Start and end of time periods 

ring ch printouts of temperature rise are desired. 
Time measured from start of exposure. For all print- 
outs set TIMEL(1)*#0, TIME2(1)=10,000 (arbitrarily 
large time), and KX=1. 
Unit: sec 

TO--Tnitial temperature of skin and subcutaneous tissues. 
Unit: °C | 


Reconmended value: 37 


.TSWEAT--Thickness of sweat layer. Sweat layer considerwd 


only if it has thickness exceeding the arbitrarily small 
value of 0.0001 em. 


Unit: om 
Recommended values: 0 to 0.01 

WAVEL--Wavelength of laser's radiation. 
Unit: pm 

WO,W1,W2,W3--Water content associated with sweat, anterior 
epidermal layer, dermis, and sut-cutaneous tissues, 
respectiveiy. 


Units: g/om? 


Recommended values: WO 


= 
NO 
HO dw H 
© 
oo 
So 


XC--Factor by which time steps are sequentially increased. 
Excessively large XC values will cause temperatures to 
oscillate. 

Unit: unitless 
Recommended values: 1.1 to 1.4 


ZBL--Depth at which steam blister forms, equais thickness of 
epidermis. 


Unit: om 


Recommended values: 0.0121 
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aDEP--Depth at which water is prevented from evaporating and 


may become superheated. 

init: oom 

Recommended value: 0.9060 (half epidermal thickness) 
ZR--Factor controlling value of smallest time step. 

Unit: unitless 

Recommended value: 0.5 


Z2W(L), =! NW--Arbitrary factors used in comouting rates of 
eat Toss due to the generation of steam. 


Unit: unitless 
Recommended value: 1 


2Z--Factor by which minimum time step is altered following 
each pulse. 


Unit: unitless 


Recommended value: 1 
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APPENDIX B 
PROGRAM LISTING 


@ Main Program 

@ Subroutine TIME 
@ Subroutine GRID 

® Subroutine PROF 

@ Subroutine CWATER 
@ Subroutine BA 

® Subroutine HTDEP 
@ Subroutine TEMP 

@ Subroutine DAMAGE 
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